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An Ha-selected sample of cataclysmic variables — II. 
Implications for CV evolution 
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ABSTRACT 

We use an independent new sample of cataclysmic variables (CVs), constructed by 
selecting objects for Ha emission, to constrain the properties of the intrinsic CV pop- 
ulation. This sample is restricted to systems that are likely to be non-magnetic and 
unevolved; it consists of 17 CVs, of which at least 10 have orbital periods above 3 h. 
We find that even very generous allowance for selection effects is not sufficient to rec- 
oncile the large ratio of short- to long-period CVs predicted by standard CV evolution 
theory with the observed sample, possibly implying that short-period systems evolve 
faster than predicted by the disrupted magnetic braking model. This would require 
that an angular momentum loss mechanism, besides gravitational radiation, acts on 
CVs with orbital periods below the period gap. To bring the model into agreement 
with observations, the rate of angular momentum loss below the period gap must be 
increased by a factor of at least 3, unless the model also over estimates the angular 
momentum loss rate of long-period CVs. 

Key words: binaries - stars: dwarf novae - novae, cataclysmic variables. 



1 INTRODUCTION 



Despite its wide relevance, the evolution of cataclysmic vari- 
, ables (CVs) is still not very well understood. Much of the 
' problem stems from the fact that strong selection effects act 
[ on observed samples of CVs, making it difficult to constrain 
1 theory observationally. 

The angular momentum loss rate (—J) is the crucial 
ingredient of CV evolution theory. Angular momentum loss 
leads to mass transfer from the secondary to the white dwarf. 
The relation between the mass transfer rate (M) and J de- 
pends on the structure of the secondary star, and the re- 
action of the secondary to mass loss determines the orbital 
period {Port) evolution of the system. 

Population synthesis methods combine the changing M 
and Porb with a model of the birthrate of CVs to predict the 
present-day distribution of CVs as a function of Mi , M2 (the 
mass of the white dwarf and secondary, respectively), M, 
and Porb- Only one of these parameters, Porb, can practically 
be m easured for a large number of CVs (but see iPattersonI 
Il998l for an indirect method to measure mass ratios). The 
orbital period distribution of known CVs is therefore one 
of the few observational properties of the CV population 
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that can be used to constrain theory. The other obvious 
constraint is the observed CV space density. 

In comparing both the size and period distribution of 
the known CV population to theoretical predictions, obser- 
vational bias must be accounted for. This requires that ob- 
served CV samples be restricted to objects selected in well- 
defined (and preferably homogeneous) ways. 

A long-standing problem is that theory predicts a large 
population of short-period CVs, consisting mostly of period 
bouncers (systems that have evolved beyond the period min- 
imum at ab out 7 6 min), which i s not observed. The models of 
iKolbl (|l993l ) and lHowell et al.l l|l997h both predict that only 
~ 1% of CVs are long-period systems, while roughly 70% are 
period bouncers. Existing observations have already been 
used to argue that the shor t-period CV population cannot 
be as l arge as predicted (e.g. |Patterson|[T998l : |Pretorius et al.l 
l2007al ). We will do the same here, using a new and indepen- 
dent CV sample. 

The CV samples that have been available to date are 
heavily biased against the intrinsically faint, short-period 
CVs, mainly because most have bright flux limits. The new 
sample considered here is also limited to apparently bright 
systems. However, almost all surveys incorporate a second 
selection cut (most commonly a blue cut) that also discrim- 
inates against the discovery of short-period CVs. The CV 
sample that we have constructed differs from most existing 
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samples in that the only selection criterion (other than a 
flux limit) is based on line emission. 

The spectra of the majority of CVs show Balmer 
emission lines, originating mainly in the accretion flow. 
The luminosity of CVs is anti-correlated to the equiva- 
lent widths (EWs) o f their emission lines (|Pattersonl 1 19841 : 
IWitham et"all l2006l : we will take EWs of emission lines 
as positive throughout). The explanation for this anti- 
correlation is that intrinsically faint CVs are low-M systems 
with low density discs, in which recombination line cooling is 
very efficient. Therefore, in contrast to other selection crite- 
ria (such as blue optical colours and variability) , an emission 
line EW-based selection cut should favour the discovery of 
intrinsically faint, short-period systems. 

We have used the AAO/UKST SuperCOSMOS Ha Sur- 
vey (SHS) to define a homogeneous sample of CVs, se- 
lected on the basis of Ha emission. Obs ervations of the new 
CVs were presented in an earlier paper l|Pretorius fc Kniggel 
l2008l : hereafter Paper I). Here we describe the construction 
and completeness of the sample, examine the observational 
biases affecting it, and compare it to the predictions of the- 
ory. 



2 IDENTIFYING CVS IN THE SHS 

The SHS is a photogra phic survey carried out with the 
UK Schmidt Telescope jParker et al.l |2005| ). Plates were 
scanned by the SuperCOSMOS digitizing machine, and the 
data are publicly available. The survey covered a total of 
^ 3 700 sq.deg. (233 separate 4° x 4° fields) at low Galac- 
tic latitude in the southern hemisphere {\b\ < -|-10° and 
5 < +2°), in R and narrow-band Ha. The limiting magni- 
tude is i? ~ 20.5, and the angular resolution of the images 
is ~ 1". 

The digitized plates are processed to extract sources 
and provide Image Analysis Mode (lAM) data for the de- 
tected objects. The analysis assigns, amongst other param- 
eters, classification, quality, and blend fiags to each source. 
This makes it possible to obtain catalogues containing only 
well-isolated point sources with relatively good photometry. 

Older UKST broad-band surveys, including an /-band 
survey of the wh ole southern sky, have also been scanned by 



survey ot tne wn oie soutnern skv, nave also Deen scannea Dy 
SuperCOSMOS ( Hamblv et al]|2001al : lHamblv et al.ll2001bl : 
iHamblv et al.. .200l3 '). The photomet r ic accu racy of these 
surveys is discussed in iHamblv et al.l (|2001bl ). The SHS R 
and Ha data are expected to be of similar quality, since 
the survey was digitized and processed in the same way. 
The photometry is calibrated to remove systematic errors 
in colour as a function of ma gnitude and position on a 
given plate; t his is described in IHamblv et al.l (l2001bl ) and 
iParker et al.] (|2005| ). The result is that errors in colours are 
typically smaller than would be expected from the errors 
in magnitudes. By matching the UKST catalogues to CCD 
data from the INT/WFC photometr ic Ha survey of the 
northern galactic plane (IPHAS; see iDrew et al.l l2005f ) in 
regions near (5 = 0° where the two surveys overlap, we made 
a rough estimate of errors in the SHS colours. Typical errors 
for stars brighter than R = 17 are ~ 0.07 in ii — Ha and 
~ 0.14 in 7? — I, although these errors vary from field to 
field. 

Given the limited quality of photographic photometry. 
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Figure 1. Colour-colour diagrams of l°xl° sections of four of the 
SHS fields in which CVs were identified. The solid lines are least- 
squares fits to all data points, while the dashed lines are vertically 
displaced by 4x the rms scatter in i? — Ha around the best fits. 
Only sources above the dashed lines are considered for selection. 
Each panel is labeled with the name of the CV selected from the 
plot, and the colours of the CVs are plotted as circled points. The 
other points above the cut-off are plate defects (crosses), objects 
that were not included in the sample, despite satisfying the cut- 
off (filled circles), and an object that was observed but has no Ha 
emission (triangle). 



we do not expect any selection procedure to identify emission 
line objects in this survey with complete reliability. Poor 
photometry is also not the only reason for the selection to 
fail. Because Ha and corresponding _R-band plates are in 
many cases taken months or even years apart, and all I- 
band images were obtained decades earlier, variable stars 
can erroneously be included in the sample. Also, since M 
stars have strong molecular absorption bands in both the R 
and / pass-bands, they are another class of objects that can 
be selected, even if no emission lines appear in their spectra. 
Multi-object spectroscopy has been used for the follow-up of 
objects from this survey (e.g. Hopewell et al. 2005), and has 
yielded low hit-rates for identifying emission line objects. We 
therefore chose selection criteria that are quite conservative, 
and that target only objects with large R — Ha excess. 



2.1 The selection procedure 

Our selection is aimed at bright objects that are clear out- 
liers in the R — Ha vs R — I colour -colour plane. The se - 
lection procedure is similar to that of IWitham et al.l (|2006D . 
but simpler (we use only a single magnitude bin, and do not 
attempt to isolate an unreddened main sequence, since this 
sequence is typically not separately discernible in the SHS 
photometry). 

We considered 175 out of the 233 fields (i.e. 2 800 sq.deg. 
in total). These fields were chosen to have relatively good 
photometry (colour-colour diagrams of 2° x 2° sections of 
fields were inspected, and in cases where one or more of these 
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sections had a very poorly defined stellar locus, the whole 
field was discarded); fields observed in the multi-object spec- 
troscopy proRram of which some results are presented by 
iHopeweU et al] (|2005l ) were also avoided. 

We extracted catalogues covering 4° x 4° around every 
of the 175 field centres, and restricted these catalogues to 
unblended point sources brighter than R = 17.0, and with 
good quality photometry in R, Ha, and jQ- 

Since the position of the stellar locus in the R — Hq 
vs R — I plane changes, even across a single 4° x 4° field 
(mainly because the effective photometric calibration is not 
constant), each field was divided into 16 1° x 1° sub-fields, 
and these 1 sq.deg. sub-fields were used individually in se- 
lecting Ha excess sources. 

Fig. [T] shows the R — Ha excess as a function oi R — I 
colour of objects in four such sub-fields. The solid lines are 
linear fits to the colours of all objects; note that the stellar 
locus has a positive slope in all these plots. The necessary, 
but not sufficient, criterion for selection is that targets have 
R — Ha > 4x the rms scatter in i? — Ha around the best 
linear fit to the stellar locus. The dashed lines are vertically 
displaced by this amount from the best fits, and are the 
cut-off criteria for inclusion in the sample. 

In fields where not all objects above the cut-off are 
cleanly separated from the stellar locus, some objects satis- 
fying the cut-off are not selected. Examples of such objects 
are plotted as filled circles in Fig. [1] This amounts to sub- 
jectively excluding some potential targets, based on inspec- 
tion of the colour-colour diagrams (or, equivalently, using a 
stricter R — Ha selection cut in some fields) . 

The colour-colour diagrams in Fig. [T] are for fields from 
which two new, and two previously known CV were selected. 
It illustrates that objects with fairly small R — Ha can be 
selected if they are blue m R — I (e.g. DV Sco), and also 
that, because of the non-zero slope of the stellar locus, a 
selection algorithm based only on R — Ha will not success- 
fully differentiate between true outliers and objects with red 
R~I. 

A selection based on abnormal colours, not surprisingly, 
identifies a large number of objects with spurious photome- 
try in one or more wave-band. Therefore, we inspected the 
Ha-, R-, and /-band images of all potential targets to ex- 
clude objects that are outliers because of, e.g., plate defects 
(some of the various ki nds of plate defects t hat affect the 
data are illustrated iu iHamblv et al. 2001a and lParker et al.l 
I2OO5I ). Outliers that were excluded from the target list for 
this reason are plotted as crosses in Fig. [1] Next, we also 
excluded blends, objects from images that appear to have 
been taken in particularly poor seeing, and a few objects 
in areas with very high Ha background. Less than 30% of 
potential targets survive these quality cuts. 

With the exception of objects that are known only as 
X-ray sources, the targets that could be matched to known 



^ To be exact, we considered objects that satisfied all of 
the following criteria (i? and Ha magnitudes are here de- 
noted SR and R_Ha, respectively). SR < 17, R_Ha<99, 
I<99, QUAL_SR<128, QUAL_Ha<128, QUALJ<128, 
BLEND.SR = 0, BLEND.Ha = 0, BLEND.! = 0, and CLASS = 2. 
The meaning of the q uality flags (e.g. QUAL_SR) is described by 
iHamblv et al] ll2001bD : the flag CLASS = 2 is assigned to objects 
classifled as stellar. 
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Figure 2. SHS colours of all objects for which we obtained identi- 
fication spectra. An estimate of typical la errors in these colours 
are shown in the top left hand corner. Several objects lie outside 
the R — I and R — Ha ranges plotted. Note that the relative posi- 
tions of objects in this plot can be misleading, since the position 
and slope of the stellar locus vary. 



objects in the SIMBAD Astronomical Database were also 
excluded. The known objects included planetary nebulae, T 
Tauri stars, CVs, and several other types of variable stars. 

The selection produced 507 targets, and we obtained 
identification spectra for 460 of them (no additional criteria 
were applied to choose these 460 objects; we simply observed 
as many targets as time allowed for). 172 of the observed tar- 
gets really are Ha emission line objects (a hit-rate of only 
37%), and 14 (that is 3%) are CVs. Of the 288 targets with 
no Ha emission, 10 are M stars (we also found many M stars 
with Ha emission). Some combination of poor survey pho- 
tometry and variability (probably in most cases continuum, 
rather than emission line variability) led to the selection of 
the remaining 278 targets found to have no Ha emission. 

Fig.[2]shows _R — Ha excess as a function of R — I colour 
for all the objects we observed. The colours of CVs are plot- 
ted as crosses; other emission line objects are shown as open 
square points, and objects without Ha emission are plot- 
ted as grey triangles. Note that some of the objects without 
Ha emission have very large R — Ha, and that the CVs are 
systematically bluer in _R — / than the overall target sample. 

We also investigated other optical colours, as well as 
near-IR colours of the targets. While it is possible to sep- 
arate the majority of late type emission line stars from 
CVs using e.g Tw o Micron All Sky Survey (2MASS; 
ISkrutskie et ahlboOSl ) photometry, Be stars and other early 
type emission line stars, as well as the objects without any 
line emission (our most important contaminant), occupy the 
same regions as most CVs in all the colour-colour and colour- 
magnitude planes. It will therefore not be easy to construct 
a much more efficient selection procedure to find CVs in this 
survey. 

In Paper I we included two CVs fainter than R — 17. 
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These were discovered during an observing run for which 
we experimented with selecting objects in a fainter magni- 
tude bin. We will here consider only CVs from the i? < 17 
selection. 



2.2 Recovery of previously known CVs 

The fraction of known CVs recovered by our selection gives 
an indication of the completeness of the CV sample. There 
are 42 previously known CVs with _R < 17 in the survey area 
covered by our follow-up, and only 7 of these were selected as 
targets. The majority of the CVs that were not selected (30 
systems) were excluded by the cuts in the class, quality, and 
blend lAM parameters. Of the other 5 that were missed, one 
was excluded because it is blended (despite being classified 
as a single point source in the SHS), and two did not satisfy 
the R — Ha ys R — I cut-off criterion. The remaining two 
CVs that were not recovered had colours above, but very 
close to the selection cut. 

With 7 recovered systems, and 14 new CVs, we know 
that our completeness is no better than (14-|-7)/(14 + 42), or 
38%. Since we observed 460 out of 507 targets, the efficiency 
of selecting previously known CVs is (7/42)(460/507). I.e., 
we expect to select roughly 15% of all CVs at i? < 17, if 
known CVs occupy the same general part of the R — Ha vs 
R — I plane as the intrinsic CV population. Fig. [3] shows the 
R, Porb, and R — I distributions of CVs with _R < 17 in the 
survey area considered, highlighting the previously known 
CVs chosen as targets by the selection. 

Although the CV recovery rate is low, a more impor- 
tant concern is magnitude- or period-dependent bias in the 
recovery rate (as that would almost certainly imply the same 
bias for the overall CV sample). It is important in this re- 
gard that a large fraction of CVs known to have been missed 
by our selection (30 out of 35) were excluded by the restric- 
tions on lAM parameters, which are expected to be mostly 
insensitive to colour and apparent magnitude. 

The magnitude distributions of both the selected CV 
sample and the recovered CVs appear to indicate that, up 
to R = 17, our survey is not affected by an apparent mag- 
nitude bias. This is illustrated both by the steep rise to- 
wards fainter magnitudes in the R distribution (left hand 
panel of Fig. |3}, and by the fact that we do not recover 
relatively more bright CVs (the same fraction of previously 
known CVs are recovered in the magnitude bins R < 14, 
14 ^ i? < 16, and 16 i? < 170 The middle panel of 
Fig. [3] also shows clearly that a larger fraction of known 
CVs are recovered at short periods than at long periods. It 
is therefore tempting to conclude that, as expected from the 
fact that short-period CVs have stronger emission lines, any 
period bias in the selection favours the inclusion of short- 
period CVs. 

The right hand panel of Fig. [3] shows the distribution 
of CVs in i? — 7. The distribution of the sample of previ- 
ously known CVs has a standard deviation of 1.0. This large 



^ Note that the magnitude distribution of the sample of pre- 
viously known CVs indicate that it is not appr oximately com- 
plete to R = 17. This is not surprising; see, e.g.. iPretorius et al.l 
ll2007al) . 



spread in R—I is not physical — it results from poor photom- 
etry and variability between the R and / epochs. The CVs 
included in our sample are on average bluer than the previ- 
ously known CVs. This bias in i? — / is not as strong as it 
appears at first sight. We will not consider CVs at very long 
periods, because they are likely to contain evolved secondary 
stars. Seven systems in the sample of previously known CVs 
have periods greater than 8 h, and 5 of these have R — I 
colours between 0.36 and 0.70, i.e. in the two maximum bins 
of the R — I distribution. With these 7 objects excluded, the 
probability that the R — I distributions of the selected and 
previously known CVs are drawn from the same population 
is 0.54, according to a Kolmogorov-Smirnov test. Therefore, 
we have no good evidence for a blue bias in i? — 7 from this 
comparison. However, the samples used for this comparison 
are both quite small, and the positive slope of the stellar 
locus in the 7? — Ha vs R — I plane, which causes the se- 
lection to be more sensitive to objects with relatively small 
EW(Ha) at bluer R — I, implies that such a bias may be 
expected. The implications that this may have for Port- or 
7?-dependent biases will be further discussed in Section [4.21 



3 THE CV SAMPLE 

The selection identified a total of 21 CVs; two thirds of this 
sample consists of newly discovered systems. Some prop- 
erties of the CVs comprising the sample are given in Ta- 
ble[T] The lower limits on the distances (d;) of the previously 
known CVs were estimated in the same way as described for 
the new CVs in Paper I. The last column of Table [T] gives 
references for the sub- types and orbital periods, where avail- 
able. 

Although we are not able to provide firm classifications 
for any of the newly discovered CVs, it is likely that the 
majority of them are NLs (see Paper I). There is no good 
evidence that any of the new systems are magnetic; however, 
they have not yet been thoroughly studied. 

V1040 Cen, CU Vel, and DV Sco are SU UMa stars. 
The orbital period listed for DV Sco is uncertain because it 
was estimated from the superhump period (the superhump 
period is 1.71 h, but even this is uncertain, because the su- 
perhump photometry available for this system is aliased). 
V1040 Cen and CU Vel both have normal fractional super- 
hump period excess, e — (Psh — Porb) /Porb, and therefore ar e 
not expected to be period bouncers ([Patterson et al.ll2003l ). 
V1040 Cen reaches V = 12.5 at maximum, and is likely to 
be very nearby. The second bright DN in the sample, V383 
Vel, was in outburst during the SHS observations. It is a 
poorly studied system; the orbital period is not known. 

The sample includes two magnetic CVs, namely VV 
Pup and HZ Pup. HZ Pup is a classical nova (nova Pup 
196 3), as well a s an IP . VV Pup is a polar. 

iKato et all l|2003D report a 37-min modulation in out- 
burst photometry of the DN CW Mon, and suggest that the 
system might be an IP. However, the high coherence associ- 
ated with the spin cycle of a white dwarf has not yet been 
demonstrated for this signal, and it is detected only near 
maximum of outburst. IWarnerl (|2004l ) classifies the 37-min 
modulation as a quasi-period oscillation — a phenomenon 
commonly seen in high-M, non-magnetic CVs. Furthermore, 
CW Mon does not have particularly strong He ii A4686 emis- 
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Figure 3. The R, Porbt a^nd R ^ I distributions of CVs with < 17 in the area covered by our follow-up. The bold histograms show 
all systems in the selected sample. Previously known systems are plotted as fine histograms, and previously known systems recovered 
by our selection are shown as cross hatched histograms. The period gap is indicated in the middle panel by light grey shading of the 
period range 2 to 3 h. One known CV in our survey area, the nova V1017 Sgr, with Porb = 5.714 d, is not shown on the scale of the Porb 
histograms. 



sion (|Szkodvll 19871 ). Note, however, that the IP HT Cam has 
a spectral appearance very similar to CW Men, and an in- 
creased white dwa rf spin pulse amplitude during outburst 
jKemp et al.ll2002h . 

Magnetic systems are not specifically dealt with in 
the population model we will use in the next sec- 
tion, and may evolv e differ entl y from non-magnetic 
CVs (e.g. IWu fc Wickramasingy Il993l : iLi et all Il994l : 
IWebbink fc Wickramasinghd 2003) ■ We will therefore ex- 
clude the magnetic CVs HZ Pup and VV Pup from the 
observed sample. CW Mon is also excluded from the sam- 
ple, despite the evidence for its magnetic nature not being 
compelling, because (in view of it being a long-period sys- 
tem) excluding it is a conservative choice. The conclusion 
regarding the ratio of long- to short-period CVs that will be 
derived in Section [5] would be even stronger had we included 
it in the sample. 

We also exclude systems with Porb > 5 h from both the 
model (see Section |4ll and observed populations, since CVs 
with evolved secondaries pr obably dominat e the populat ion 
at these periods tBeuermann et al ] 1 19981 : Haraffc & Kolbl 
I2OO0I : iPodsiadlowski et al.ll2003h . 

Therefore VV Pup, HZ Pup, CW Mon, and Ha074208 
are not considered further. The entries for these 4 systems 
are in italics in Table [1] The rest of the sample (17 CVs, of 
which at least 10 are long-period systems) likely consists of 
non-magnetic, unevolved CVs. There is no indication that 
any of these CVs are period bouncers. 



4 MODELLING SELECTION EFFECTS 

We use a Monte Carlo technique to quantify the effects of 
the most important selection biases acting on our sample. A 
random sample of CVs is drawn from a predicted intrinsic 
CV population, and distributed in a model Galaxy. By mod- 
elling the spectral energy distribution (SED) and outburst 
properties of each CV, we predict the observational appear- 
ance of this intrinsic population in a survey with our par- 
ticular selection criteria. We use the computational method 



of IPretorius et ahl (|2007al ). specifically, their model Al, with 
only a minor change to the SED model. 

The Monte Carlo method is described in Section 14.1 
below. It does not account for all the selection biases present 
in our sample. In Section 14.21 we discuss the assumptions 
that were made to deal with the remaining biases. 



4.1 The Monte Carlo code 

The Monte Carlo calculation uses as input a probability dis- 
tribution function (PDF) for the present-day population of 
CVs over the parameters Porb, {M) (the secular average M), 
Ml, and M2. Thi s PDF results from the population synthe- 
sis model p m5 of iKolbl (119931. It relies on the predicted CV 
birth rate of de Koo J 1 1992) an d the magnetic braking pre- 
scription of IVerbunt fc Zwaaiil (|l98ll ). 



The Galaxy is treated eis an axisy mmetric disc with 
a Ga lactic Centre distance of 7620pc (jEisenhauer et ahl 
I2OO5I ), a radial scale length of 3 000 pc, and vertical scale- 
heights of 120, 260, and 450 pc for long-period CVs, pre- 
period minimum short-period CVs, and period bouncers, 
respectively. The assumption of different scaleheights for 
these three sub-populations of CVs is motivated by their 
different typical ages, and are chosen to be representative 
of stellar populations with age s of 10^ y. ~ 10 ^'^ y, 
and ~ 10^ '' y, resp ectively (|Robin fc Crezel 1 19861 '): see 
IPretorius et al.l l|2007ah for a more complete discussion, and 
for an indication of the sensitivity of the period distribution 
of magnitude-limited samples to this assumption. 

We consider fiux from the accretion disc, bright spot, 
white dwarf, and secondary star in the model of the over- 
all SED of each CV in our simulation, wit h all except the 
secondary treated in the same way as in IPretorius et al.l 
(2007^). The accretion disc contribution is obtained from 
the model of iTvlendal { 198ll ). the bright spot is modelled 
as a black body, and the white dwarf fiux is found from 
the tabulation of lBergeron et al.l lll995l). using the effective 
temperature predicted bv lTownslev fc Bildsteii! (|2003l ). The 
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Table 1. Sub-types, orbital periods, _R-band magnitudes, and 
lower limits on distances of the CVs included in the survey. En- 
tries for the systems that will not be included in the comparison 
to theory are in italics. 



Object 


Type 


Porb/h 


R 


cZ;/pC 


References 


V1040 Cen 


su 


1.45 


13.9 


40 


1 


DV Sco 


su 


1.65: 


16.7 


270 


2,3 


VV Pup 


polar 


1.67392 


14.6 


140 


4,5 


CU Vel 


SU 


1.88; 


15.8 


120 


6 


Ha092134 




3.043 


16.8 


530 


7 


Ha073418 


SW: 


3.18542 


16.7 


320 


7 


Ha074655 




3.3984 


14.0 


180 


7 


Hal02442 




3.673 


16.4 


580 


7 


Hall2921 


NL: 


3.6851: 


15.5 


510 


7 


Hal03135 


NL: 


3.76: 


16.3 


630 


7 


Hal03959 


NL: 


3.785: 


16.4 


670 


7 


Hal30559 




3.93: 


16.5 


500 


7 


Ha092751 


NL: 


4.1 


16.2 


770 


7 


CW Mon 


DN,IP: 


4.238 


15.4 


210 


8 


Ha094409 


SW: 


4.506: 


16.2 


520 


7 


HZ Pup 


CN,IP 


5.11 


16.5 


1010 


9 


Ha074208 




5.706 


15.3 


490 


7 


Ha075648 


SW: 




16.0 




7 


Hal63447 


NL: 




16.7 




7 


Hal90039 


NL: 




16.8 




7 


V383 Vel 


DN 




12.3 
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References: (11 [Patterson et al.l ||2003|) : (2) v snet-alert 8325 ; (3) 
Daisaku Nogami, pr i vate c omm unication: (4) [Herbi ej ll96d|); (5) 
[Schneider fc Younel lll980l); (6) iMennickent fc Diaj (Il996l): (7 ) 



Pap er I: (8) i Kato et al.l ll2003h : (9) lAbbott fc Shafted lll997^ ; 
(10 IWilliamsl ||2000^ . 

Notes: The types are dwarf nova (DN); nova like variable (NL); 
SU UMa star (SU); SW Sex star (SW); intermediate polar 
(IP); classical nova (CN). Uncertain values or classifications are 
denoted by ':'. Hal90039 is classified as a CV based only on a 
single, low quality spectrum. 

flux contribution of the secondar y star is obtain ed from the 
semi-empirical donor sequence of lKnigg3 l|2006l 'FI. 

The only type of variability included in the model is 
DN outbursts. The modelling of outbursts relies on the 
well known empirical relation between Pprb and the a bso- 
lute magnitude of DNe at maximum (e.g. lWarneJll987l ): see 
IPretorius et~al] (|2007al ) for a more complete discussion. The 
probability of finding a DN in outburst in a single-epoch 
survey is then given by the outburst duty cycle (C). 

The R and / magnitudes produced by our SED model 
are on the Kron-Cousins RI system, and were trans- 
formed to the photogr aphic system of t he survey using 
[Morgan fc Parked (|2005l ) and lBesselll l|l986l ). The differences 
between the photometric systems are small. 

We then select model CVs with R<17 and |&| < 10 deg. 
Fig. [4] shows the differential and cumulative apparent mag- 
nitude distribution, as well as the Porb distribution, of the 
CV population predicted to be observed in our survey, if we 



^ We interpolated the high-resolution version of the donor se- 
quence onto the period grid of the model CV population, where 
possible. The upper edge of the model period gap is slightly above 
the observed location. The model period gap is also narrower than 
observed, so that we had to extrapolate the empirical sequence 
at the lower edge of the gap. 



neglect other selection cuts. The grey histogram in the left 
hand panel is the R distribution of the 17 observed CVs. 
The middle panel shows the cumulative R histograms of 
the model (black) and real samples (grey). A Kolmogorov- 
Smirnov test gives the probability of the model and observed 
R distributions being drawn from the same parent popula- 
tion as 0.35, implying that the magnitude distributions are 
consistent. In the right hand panel, the bold histogram is 
the Porb distribution of all systems detected; the contribu- 
tion of period bouncers is shaded in dark grey, and the fine 
black histogram shows the distribution of outbursting DNe 
expected in a single-epoch survey. The model histograms are 
scaled as described in Section [5.21 below. 



4.2 Additional biases 

The above modelling accounts for the effects of the magni- 
tude limit, and the restriction of the survey to the Galactic 
plane. With only those selection effects included, the pre- 
dicted period distribution is inconsistent with the observa- 
tions (relatively too many short-period CVs are predicted). 
However, we need to determine whether the remaining se- 
lection effects present in our sample could reconcile the pre- 
dicted and observed distributions. We therefore make con- 
servative assumptions about the infiuence of these biases 
(i.e. assumptions that decrease the predicted ratio of short- 
to long-period CVs). 



4.2.1 The 'EiWiYLa)- dependent selection cut 

The dependence of selection probability on EW(Ha) is not 
included in our model. This is because the disc model pre- 
dicts EWs that arc clearly too small (as demonstrated for 
H/3 by. Patterson 1984: this is also true for Ha). The alterna- 
tive of using an empirical relation is not available, because 
the dependence of EW(Ha) on, e.g., My is not known. One 
may expect that including this selection cut in the model 
would increase the predicted ratio of short- to long-period 
CVs, and therefore that not accounting for it is a conserva- 
tive assumption. 

Whether an object is selected as a target depends on its 
J? — Ha and R — I colours. The R — I dependence arises from 
the non-zero slope of the stellar locus in the J? — Ha vs R — I 
plane. It is also likely that the selection probability depends 
on R magnitude, since the photometric uncertainty increases 
at fainter magnitudes. While our CV sample does not at 
first sight appear to be increasingly incomplete at fainter 
R, the sample of all emission line objects from our follow- 
up certainly is magnitude biased, with an R distribution 
that flattens off at around R — 16. The reason that the 
CV sample is not affected by the same magnitude bias as 
the overall emission line sample is not because CVs have 
systematically stronger emission lines than other types of 
emission line objects selected for follow-up, but because the 
CVs are on average bluer in i? — /. Our selection is sensitive 
to faint objects with small EWs only if those objects are 
sufficiently blue in R — I. Redder faint objects must have 
stronger lines in order to be selected. 

To the extent that data can guide us, the result- 
ing blue bias does not matter — intrinsically fainter, red- 
der CVs have stronger lines. This is well established for 
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Figure 4. The differential and cumulative magnitude distributions, as well as the differential orbital period distribution, predicted by 
the simulation. The bold histograms in all three panels are the distributions of all systems predicted to be detected, with only the flux 
limit and \b\ constraint imposed. The R and cumulative R distributions of the 17 observed systems are also shown, in light grey, in the 
left hand and middle panels. In the right hand panel, the contribution of period bouncers to the orbital period distribution is shaded 
in dark grey, and the fine histogram shows DNc that are found in outburst in a single epoch. The model R. and Port distributions are 
scaled as described in Section 15.21 



normal short-period CVs, but observational information on 
period bouncers is scarce. There are now two confirmed 
period bouncers, SPS S J103533.03+055158.4, and SPSS 
J150 722. 30+523039.8 (jLittlefair et al.ll2006l : iLittlefair et all 
I2OO7I : Patterson, Thorstensen & Knigge (in preparation)) 
and they have i5W(Ha) = 44 A and 78 A, and Kron-Cousins 
- 7 of 0.03 and 0.04, respectively ijSzkodv et all l2005l: 
Szkodv et al ] |2006l : colour transformations fromlJester et al.l 



20051 ). Our selection is easily capable of identifying objects 



with these EWs and colours near R — 17. 

In the absence of an undetected population of faint, red 
CVs with weak emission lines, our EW cut thus preferen- 
tially excludes intrinsically bright (i.e. long-period) CVs. 



4.2.2 Variability 

PNe in outburst probably ha ve smaller EW( Hq) as a rule 
(this is certainly true for H/3: |Patterson|[l984l '). This means 
that an outbursting PN is less likely to be identified in an 
Ha survey than a quiescent system. The problem here is 
slightly more complicated than whether a system has suffi- 
ciently large EW(HQf) at a single epoch, since observations 
in the three wavebands used in the selection were not taken 
simultaneously. Variability between the epochs of the R, I, 
and Hof observations can lead both to an object with Hq 
emission not being selected, and to an object without Hq 
emission being included. The most common large amplitude 
variability displayed by CVs is PN outbursts, and we will 
consider only this type of variability. Out of the 175 fields 
we used in the selection, 72 were observed in R and Hq on 
the same night; all /-band data were taken roughly 20 years 
earlier. 

A system that is bright in both R and Hq but faint in 
/ will probably be selected, because it will be very blue in 
R~ I (V383 Vel is an example — it was observed in outburst 
in R and Ha on the same night, but was faint during the I 
epoch). CVs are likely to be missed by the selection because 
of variability if they are bright in R and faint in Ha, or 
bright only in /. We assume that the probability of a long- 



period system being excluded for reasons of variability is 
0, and that the probability of a short-period system being 
excluded for being bright in R, is C. Then, to allow for the 
possibility of a PN being bright in /, we do the same with 
the remaining short-period PNe (i.e. we exclude a fraction 
of systems equal to the particular outburst duty cycle). 



5 COMPARISON TO THEORY 

5.1 Relative numbers of short- and long-period 
CVs 

After applying the cuts involving R, \b\, and variability de- 
scribed above, the sample is predicted to contain 13% period 
bouncers (i.e. about 2 systems in our sample of 17), whereas 
none are detected. This is not a serious discrepancy. The 
predicted fraction of long-period CVs is 30%. If we assume 
that our combined EW and R — I cut makes it impossible to 
detect period bouncers, then 34% of the sample is predicted 
to be long-period systems. Using the binomial distribution, 
the probability that the observed ratio of at least 10/17 is 
consistent with this prediction is 0.031. The model is there- 
fore inconsistent with observations at more than 2-a. 

The inconsistency would be a 3-cr result if two more 
systems in our sample had periods above the period gap. It 
is likely that this will prove to be the case. The spectrum 
of Hq190039 is characteristic of an intrinsically bright CV, 
and there is reason to suspect that Hq075648 is an SW Sex 
star. 



5.2 Space density 

Our sample is not suitable for measuring the CV space den- 
sity (p), but we can verify that it is i n reas onable agreement 
with a specific prediction. Ide Kooll l|l992t ) a ndlKolb (19931) 
theoretically predict p ~ 10~* (but see also IPolitanol [l996l . 
for a theoretical value significantly lower), with a hundredth 
of all CVs above the period gap. The space density of long- 
period CVs (pi) is then predicted to be ~ 10~®pc~^. 
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Setting pi — 10~® pc~'^ in the simulation, we predict 
(with the area coverage and flux limit of our survey) that 
about 22 long-period CVs should be detected if the sample 
really is 15% complete. This assumes that the only selection 
cut affecting long-period CVs is the magnitude limit and 
|6| restriction. Our sample contains at least 10 long-period 
CVs, perhaps as many as 14. The total number of observed 
long-period CVs is then consistent with the prediction if the 
true completeness for long-period systems is ~ 8%. This is 
entirely reasonable, but given the uncertain completeness, it 
is not a very secure resuHQ. 

The predicted R and Port histograms in Fig. |4] were 
scaled to reflect the number of CVs that should have been 
detected in our survey (2 800 sq.deg., |&| < 10°, and R < 17) 
if pi = 10~® pc~^, and if the survey is 8% complete, with no 
dependence of completeness on Port- 



5.3 Evolutionary time-scales 

It is usually assumed that the theoretical ratio of the num- 
ber of short- to long-period CVs is larger than observed be- 
cause theory predicts too large a population of short-period 
CVs. Several solutions to this problem have been proposed. 
We will consider here the suggestion that there is an angu- 
lar momentum loss mechanism that increases the absolute 
value of J of short-period CVs above that of the r ate result- 
ing from gravitational radiation (Jgr)', see e.g. Patterson! 
1 1998') : Patt erson (20 01); King ct al. (2002); Patterson e t ail 
l 2003l ): lBarker fc Kolbl ((200^). Higher -j below the period 



gap decreases the evolutionary time-scale of short-period 
CVs and therefore reduces the total number of CVs in this 
evolutionary phas^fl. 

The evolutionary time-scale of a CV is 



J_ 



\ M I 



where the angle brackets denote the secular average. Assum- 
ing that all CVs form at periods longer than ~ 5 h, and only 
evolve through the period ranges that we consider, the num- 
ber of CVs that exist in some particular period interval is 
proportional to the time CVs spend in that interval. The ra- 
tio of the numbers of long- and short-period CVs is therefore 
proportional to the ratio of their evolutionary times-scales 
(ji and Ts for long- and short-period CVs, respectively). 

The ratio of long- to short-period CVs observed in a 
survey for which the most important bias is a flux limit, can 
be written as 



(^) -K-r^IrM^ 
\nsJ obs Vr. 



(1) 



Here, > 1 is a monotonically decreasing function of ti/ts, 
since apparent brightness depends on accretion luminosity, 
which in turn scales inversely with r. The form of k means 



^ The available observational estimates of pi are ~ 2. 5 X 10 ^ 
a nd ~ 8 X l O "'^ l|Pattersonlll99dlPretorius et al.ll2007bh . 
* iPattersonI lll998l) also points out that a population synthesis 
model using such a higher — J could predict the position of the 
period minimum correctly, and would fit the observed relation be- 
tween e and Port better than models assuming only gravitational 
radiation. 



that, by comparing the observed and predicted number ra- 
tios, we can set a limit on the true ratio of evolutionary 
time-scales of CVs above and below the period gap. Using 
equation [T] we can write this ratio as 



pred 



k{{ri/Ts)pred) 

k{ri/Ts) 



{ni/us 



ini/ns)p^^a 



(2) 



where {ti/ts)^^^^ is the predicted ratio of the evolutionary 
time-scale of long-period CVs to the evolutionary time-scale 
short-period CVs. If (ni/na)^^^ > {ni /ns)^,^^^, as in our case, 
we must clearly have ti/ts > (''"i/''"s)pred- monotonic 
nature of k then implies that k{{Ti/Ts)pred) > k{Ti/Ts), so 
that 



HyflL) 

Ts \Ts / 



Ts / pred 



(ni/ris)^,, 



(3) 



Physically, this inequality holds because any adjustment to 
the predicted evolutionary time-scales that causes the short- 
period phase to be shorter relative to the long-period phase 
will also make short-period CVs brighter relative to long- 
period systems (and therefore less under-represented in a 
flux- limited sampl^l). The impact of such a time-scale ad- 
justment on the observed number ratio is therefore always 
smaller than on the intrinsic number ratio. 

To be more specific, we can assume iji)^^^^ = ti. Then, 



for this survey, with (rii/ns) 



pred 



3/7 and {ni/us 



10/7, the model evolutionary time-scale of short-period CVs 
must be decreased by > 3 to match the observations. Or, 
equivalently, we require — J > —3Jgr for short-period CVs, 
in order to bring theory and observations into agreement. 
This factor of at least 3 would be an under estimate if the 
EW-based selection cut discriminates strongly against long- 
period CVs. 



6 DISCUSSION 

Despite the very conservative assumptions outlined in Sec- 
tion 14.21 the model CV population and the observed sam- 
ple do not agr e e. We have thus confirmed the result of 
IPretorius etM l|2007ah with an independent observational 
sample. Standard CV evolution theory predicts too large a 
ratio of short- to long-period CVs. This is true even if we 
assume that period bouncers are undetectable. 

The favoured explanation for the large predicted ra- 
tio of short- to long-period CVs is that short-period sys- 
tems evolve faster than predicted by the disrupted magnetic 
braking model, so that the predicted population of short- 
period CVs is too larg e (e.g. iPattersoiill 19981 : |Pattersonll200ll : 
[Patterson et al.ll2003l ). Using the Ha sample, we find that 
the evolutionary time-scale of short-period CVs in our model 
should be decreased by a factor of at least 3 in order for the 
model to agree with observations (if the predicted evolution- 
ary time-scale of long-period CVs is correct) . 

For the Palomar-Grce n (PC) Survey CV sample 
(jCreen et al.| [r986: Ririgwald 19931 ) . IPretorius et all (|2007ah 
found (ni/us) — 7/93, compared to an observed ratio 



^ The EW cut will not select against short-period CVs, unless 
their luminosity becomes similar to that of long-period CVs — 
this means an increase by > 10. 
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of {ni/na)^^^^ 14/13. With the same reasoning as in Sec- 
tion [5]3] above, this imphes that the evolutionary time-scale 
of short-period CVs in the model is too long by more than 
an order of magnitude. This is surprisingly large , especially 
when compared to the result of IPattersonI (1200 ll ). who finds 
that J = 3Jgr for short-period CVs provides a satisfac- 
tory fit to his data. This may be a hint that the assump- 
tion (Ti)p^^^ = n is wrong. The observational constraints on 
the size of the long-period CV population and M above the 
period gap are not sufficiently precise to answer this ques- 
tion definitively, but it ce rtainly has been suggested that the 
IVerbunt fc ZwaanI (|l98ll ) magnetic brakin g law used in the 
model is too strong (e.g. Ilvanova fc Taami r2003'l. Our model 
may therefore not only predict too many short-period CVs, 
but also too few long-period CVs. 

In this regard, it is also worth noting that recent esti- 
mates of M in a small sample of short-period CVs are consis- 
tent with the values expected fro m angular momentum loss 



through gravitational radiation (iLittlefair et al 



l2006l : Lit- 
200'il ) have 



tlefair, et al. (in preparation')'). ILittlefair et al.l 
also estimated M for SDSS J150722. 30-^523039. 8, but they 
find that this system probably has an unusual evolution- 
ary history (see also Patterson, Thorstensen & Knigge (in 
preparation)). 

Before an Ha-selected CV sample can be used to place 
tighter constraints on CV evolution, EW(Hq) must be mod- 
elled in terms of the parameters describing CVs in a theo- 
retical population. It is known that the EWs of disc emis- 
sion lines increase with decreasing luminosity. However, CV 
discs are observed together with a white dwarf primary. Very 
faint systems are white dwarf dominated (i.e., the continuum 
emission from the disc is low relative to that from the white 
dwarf). Also, the disc emission lines are superimposed on 
the absorption lines of the white dwarf. Therefore, at faint 
luminosity, th e relation betwe en EW and My is expected 
to flatten off (|Pattersonl [iQsi plots E'W(H/3) vs the abso- 
lute magnitude of the disc alone; his relation therefore does 
not show this effect). This is observationally clear from the 
spectra of the f aintest known CVs, a nd ha s already been 
remarked on by lAungweroiwit et all l|2006l '). Although the 
decrease in EW(Hq) is expected to be less severe than that 
in EW(H/3), this might ultimately place a limit on what can 
be achieved by an emission line survey. 

The comparison of this Ha-selected CV sample to the- 
ory is also complicated by the possibility of variability be- 
tween non-simultaneous observations in the three wave- 
bands used for the selection. Surveys such as IPHAS will 
overcome this difficulty. IPHAS will also provide much bet- 
ter photometric accuracy than can be achieved photograph- 
ically, so that the photometry will allow for the identifica- 
tion of objects with smaller EWs in all parts of the colour- 
colour plane; any blue bias will therefore be weaker (see 
Drew et al. ll2005l:IWitham et al. I I2OO6I : IWitham" et al. I I2OO7I : 
Witham et al mm . Furthermore, the emission line samples 
constructed from IPHAS data will be much deeper. This is 
very important, since a selection technique capable of iden- 
tifying intrinsically faint CVs is of limited value in a survey 
with a bright flux limit. The main challenge for IPHAS will 
be the amount of observational effort needed to flnish the 
identification and follow-up of the CV sample. The same is 



true for other large, deep surveys — the completion of the 
resulting new CV samples are still some time in the future. 

With the smaller CV samples available at the moment, 
it is already possible to recognize the qualitative changes 
that are needed to reconcile theory with observations. 



7 CONCLUSIONS 

We have compared a homogeneous CV sample, selected for 
Ha emission, to a model CV population based on standard 
CV evolution theory. The magnitude limit and Galactic lat- 
itude range of the observed sample was modelled in some 
detail, while conservative assumptions were made to account 
for the effects of variability and the EW(Ha)-based selection 
cut. The model population is inconsistent with the observed 
sample. Specifically, the model predicts relatively too many 
short-period CVs. This confirms earlier results, based on in- 
dependent observations. 

The reason for the mismatch between the predicted and 
observed ratio of short- to long-period CVs may be that 
the theoretical evolutionary time-scale for CVs below the 
period gap is too long. A (very simplistic) consideration of 
the relative numbers of long- and short-period CVs included 
in the sample indicates that the disrupted magnetic braking 
model underestimates — J of short-period CVs by a factor of 
at least 3, assuming that the model is correct for long-period 
CVs. 

Although surveys now in progress will in future provide 
much better observational constraints on CV evolution the- 
ory than can be derived at the moment, it is already clear 
that the standard magnetic braking model is in need of re- 
vision. Furthermore, it seems that the correct approach to 
take is to investigate angular momentum loss rates in excess 
of the gravitational radiation rate in CVs below the period 
gap. 
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